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Context. A noticeable fraction of subdwarf B stars shows either short-period (p-mode) or long-period (g-mode) luminosity variations, 
with two objects so far known to exhibit hybrid behaviour, i.e. showing both types of modes at the same time. The pulsating subdwarf 
B star V391 Pegasi (or HS 2201+2610), which is close to the two known hybrid pulsators in the logg-Teff plane, has recently been 
discovered to host a planetary companion. 
, Aims. In order to learn more about the planetary companion and its possible influence on the evolution of its host star (subdwarf B 

^ |. star formation is still not well understood), an accurate characterisation of the host star is required. As part of an ongoing effort to 

I ' significantly improve the asteroseismic characterisation of the host star, we investigate the low-frequency behaviour of HS 2201+2610. 

O ' Methods. We obtained rapid high signal-to-noise photometric CCD (B-filter) and PMT (clear-filter) data at 2 m-class telescopes and 

carried out a careful frequency analysis of the light curves. 
t/^ . Results. In addition to the previously known short-period luminosity variations in the range 342 s - 367 s, we find a long-period 

^ , ' variation with a period of 54 min and an amplitude of 0. 15 per cent. This can most plausibly be identified with a g-mode pulsation, so 

that HS 2201+2610 is a new addition to the short list of hybrid sdB pulsators. 

Conclusions. Along with the previously known pulsating subdwarf B stars HS 0702+6043 and Balloon 090100001 showing hybrid 
^ ■ behaviour, the new hybrid HS 2201+2610 is the third member of this class. This important property of HS 2201+2610 can lead to a 

^.f^ • better characterisation of this planet-hosting star, helping the characterisation of its planetary companion as well. Current pulsation 

' models cannot yet reproduce hybrid sdBV stars particularly well and improved pulsation models for this object have to include the 

j^-^ , hybrid behaviour. 

, Key words, subdwarfs - Stars: horizontal-branch - Stars: oscillations - Stars: individual: HS 2201+2610 

o 

1 . Introduction acoustic pressure waves are of low degree and low radial order. 

^— ^ ' In contrast, the g-mode pulsators have even lower amplitudes 

^ • Subdwarf B stars (sdBs) are evolved objects with masses around (f^^ ^^^ ^ -^^^ ^30 90 min) a t lower temper- 

. ^ , half a solar mass and are thought to be core helium burning. ^^^^^^ ^^j^^ ^1^^^ prototyp es are EC 14 026-2647 (iKilkennv et al.l 

. Having lost most of their original hydrogen envelope in earlier 1997) and PG 1716+426 (G reen et al.ll200? ). respectively. Both 

H stages of their evolution, the remaining H-rich envelope is too ^^^^ ^^-^^^ ^ t he same /.-mechanism 

. thin to sustain hydrogen shell burning. The sdBs populate the due to an opacity bump of iron and nickel (ICharpinet et al.ll997t 

extreme horizontal branch (EHB) at effective temperatures of [Fontaine et al.l200llJefferv & SaiaiOOeh . Hybrid sdB pulsators 

20 000 to 40 000 K and surface gravities log(g/cm s ^) between ^j^^^ ^^^^^^ simultaneously and are of particular inter- 

5.0 and 6.2. Instead of evolving towards the asymptotic giant ^-^^^ jj^^ ^^^^ jyp^^ p^^^^ ^^^^^^^ ^^^-^^^ ^-^^^^ jj^^ 

branch, the sdBs are predicted to follow tracks leading direcdy j„ ^j^^ ^ _ j-^^ diagram, the hybrids are located at the in- 

towards the white dwarf region after leaving the extreme hori- ^^^^^^ jj^^ ^^^^^ g- mode instability regions. The known 

zontal branch. A fraction of the sdB stars shows pulsations: the ^^^^^-^^ p^S 0702+6043 (Schuhet al."2()05';2006'; 'Lutz et al. 

first pulsating sdB star was discovered in 1997. Variable subd- 20O8) and Balloon 090 100001 (Oreiro et al. 2005; Bar an et al 

warf B stars (sdB V stars) can be divided into the classes of rapid |2005|) 
p-mode pulsators (sdBVr) and slow g-mode pulsators (sdBVs), 

with two objects known so far to belong to both classes simulta- ™s P^Per presents the third member of this group: 

neously (hybrid pulsators, sdBVrs). 2201+2610 is a 14.3 B-magnitude star having an eff^ective 

The p-mode pulsators show low amphtudes (few ten mmag) temperature of 29 300±500 K and a surface gravity log g of 

and short periods (few minutes) at higher temperatures. The 5.40+0.10 m cgs umts (| 0stensen et al.| | 2001|) . Based on these 

spectroscopic parameters, which are very similar to those of 

* Based on observations collected at the Centro Astronomico HS 0702+6043 and Balloon 090100001, HS 2201+2610 was a 

Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- natural further target to search for simultaneous g-modes in a 

Planck Institut fur Astronomic and the Instituto de Astrofi'sica de known p-mode pulsator after the discovery of the first hybrids. 

Andalucfa (CSIC). The high signal-to-noise observations initiated to check for hy- 
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Fig. 1. Light curves from the September/November 2006 runs at Fig. 2. Lomb Scargle Periodograms (LSPs) of the p-mode 
Calar Alto. The relative flux is set to 1.0 for the 27. Sep 2006 regime. The top panel shows the complete Calar Alto data. The 



light curve and previous (following) nights are shifted upwards 
(downwards) by 0.1 units. 



bottom panel displays the LSP after simultaneous subtraction of 
five frequencies. Horizontal dotted lines indicate 4 cr confidence 
levels and the marks correspond to the IDs in Table [T] 



brid behaviour in HS 2201+2610 fitted in well with our ongoing 
long-term monitoring in order to search for secular changes in 
the pulsatio n periods i ndicative of slow evolutionary changes 
in the star jSilvotti et al. 2002), which has recently resulted 
in the surprising detection of a planetary compa nion around 
HS 22 01-1-2610 (better known as V391 Pegasi b. .Silvotti et all 
l2007h . In fact, from a few of the better-quahty data sets obtained 
within this overall observational effort up to the end of 2005, we 
already suspected the presence of additional power in the low- 
frequency domain of HS 2201+2610. 

In Section|2]we present the photometric data obtained specif- 
ically to investigate the gravity mode frequency range, as well 
as suitable confirmation data obtained as part of our long-term 
monitoring. The data reduction and frequency analysis are de- 
scribed in Section [3] and as a result of our investigation we 
will present HS 2201+2610 as a new member of the group 
of hybrid pulsating sdB stars. We discuss the implications of 
this discovery for the ongoing asteroseismic characterisation of 
HS 220 1 +26 1 0, in the context of our efforts to constrain the mass 
of its planetary companion, in Section |4] Finally, we present a 
concluding outlook in Section|5] 

2. Observations 

The light curves (see Fig.[T]) were obtained at Calar Alto 
Observatory in September and November 2006. The instrumen- 
tal setup was chosen to consist of the Calar Alto Faint Object 
Spectrograph (CAFOS) focal reducer attached to the 2.2 m tele- 
scope. A 2k X 2k CCD chip recorded the science frames through 
a Johnson B filter with a 2 x 2 binning and an exposure time of 
ten seconds with an average cycle time of 25.5 seconds. In order 
to sample the short-period /5-modes as well as to detect very low- 
amplitude ^-modes, one has to find a good compromise with the 
exposure time. Ten seconds turned out to be sufficient to resolve 
the p-mode cycles as well as provide a S/N good enough to re- 
solve low-amplitude signals. We have a total of about 26 hours 
of photometric data. 

The data reduction was carried out usin g the Time Resolve d 
Imaging Photometry Package (TRIPP, see ISchuh et alj |2003[) . 
which is designed to handle the aperture photometry of large 



CCD sets. This includes bias-, flatfield- and background correc- 
tions and finally yields extracted light curves in the form of rela- 
tive flux (in comparison to four stable reference stars) as a func- 
tion of geocentric Julian date. Additionally, a barycentric time 
correction was applied. Differential extinction effects on a night- 
by-night basis were corrected using second-order-polynomials. 

In addition to our Calar Alto data from 2006, one of us (R. S.) 
provided older data from August/September 2003, which were 
part of the Whole Earth Telescope extended coverage 23 (WET 
Xcov23). The data taken at Loiano Observatory (LOI, 1.5 m 
telescope and 3 channel photometer) and Observatoire de Haute- 
Provence (OHP, 1.9 m telescope and 3 channel photometer) were 
provided as extracted light curves (total length of 28 hours) and 
we used them to cross-check our Calar Alto results. 

3. Results 

We carried out a frequency analysis on both data sets with a pro- 
gram called sinf it, which is part of TRIPP. This includes least- 
squares non-linear sinusoidal fits to the light curves and calcu- 
lation of the corresponding periods, amplitudes and phases. The 
amplitudes are given in milli modulation intensity units (mmi). 
In this context 10 mmi means that the luminosity has changed 
by one per cent with respect to the mean brightness level. We 
used the Scargle algorithm ( Scargle 1982) to calculate Lomb- 
Scargle-Periodograms (LSPs) in the form of power spectral den- 
sity (psd) as a function of frequency, together with confidence 
levels, whic h are based on false a larm probabilities (faps). See 
Section 4 in iDreizler et all (12002 ) for a more detailed descrip- 
tion. In our case, we derived the faps from 158000 white noise 
simulations. This yields a confidence level of 4 cr, judging a sig- 
nal above this level to be real with a probability of 99.99366%. 

Table [T] summarizes the results of our frequency analyses of 
both data sets. We can clearly recover the three strongest m odes 
found in an extensive run published by lSilvotti et alJ (l2002h . here 
labeled with fl to f3. These frequencies are assigned to pressure 
driven modes. The original and residual LSPs of our Calar Alto 
data in the /?-mode regime are displayed in Fig.|2] In addition, we 
find a long-periodic luminosity variation (f7, shown in Fig. [3]) of 
about 54 minutes. 
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Fig. 3. Results of our frequency analysis in the g-mode regime. Left: September/November 2006 Calar Alto data. The top panel is a 
Lomb Scargle Periodogram of the three best nights of the full light curve where five frequencies have already been subtracted (f 1 to 
f5, see Table[Tli. Fitting a long-periodic luminosity variation (f7) to this subset of the data results in the residual LSP displayed in the 
middle. The window function is shown in the bottom panel. Right: August/September 2003 OHP/LOI data. The top panel displays 
the Lomb Scargle Periodogram of the light curve after pre-whitening of three frequencies (fl to f3, see Table [TJ. The residual LSP 
after subtraction of the long-period variation f7 is shown in the middle and the window function is illustrated in the bottom panel. 
Horizontal dashed lines indicate 4 cr confidence levels. 



As mentioned before, we looked at suitable WET Xcov23 
data with the focus on the g-mode domain to confirm our detec- 
tion of a long-periodic luminosity variation found in the Calar 
Alto data (left plots in Fig.|3ll. The periodogram of these addi- 
tional data is shown in the right plots of Fig.|3] The OHP/LOI 
results (also in Table (TJ clearly confirm the detection of a low 
frequency oscillation in HS 2201-1-2610. 

As can be seen in the middle right plot in Fig.|3] there is 
still further residual power above 4 cr in the OHP/LOI data LSP 
after subtraction of the long -period variation fl. However, this 
is not seen in the Calar Alto data (which were taken three years 
later) so that further high-quality data are required to test the 
persistence of these residual peaks. 

4. Discussion 

We identify the long-period low-amplitude luminosity variation 
in HS 2201-1-2610 with a gravity driven mode. The presence 
of both types of modes within HS 2201-1-2610 would therefore 
define this object to be a member of the hybrid sdB pulsator 



group, beside the previously known hybrids HS 0702+6043 and 
Balloon 090100001. Other interpretations for the long-period 
variation can be excluded with a hig h probability wit h simi- 
lar arguments as for HS 0702+6043 (ISchuh et all 12006). The 
same long-period variation, in frequency as well as amplitude, 
is present in data sets three years apart. This rules out sky or 
transparency variations, as well as rotating star spots. The fast 
rotation rate that would be required for the latter interpretation 
is further unlikely since rotation periods of single sdB stars are 
typically much longer. Rotational splitting (for / > modes) 
has not yet been detected in HS 2201+2610, but the split com- 
ponents often have lower amplitudes and might be below our 
detection limit. With the light-travel time effect used in the de- 
tection of the planet candidate around HS 220 1 +26 1 in its wide 
orbit, we would not be sensitive to a further hypothetical com- 
panion in a very close orbit with a 54 min period. A binary ori- 
gin of the long-period variation is unlikely for the following rea- 
sons. A 54 min orbital period would be at the extreme low limit 
of the observ ed and theoretical orbital pe riod distributions (see 
e.g. Fig. 2 in iMorales-Rueda et al.ll2006l) . Additionally, a low- 
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Table 1. Results of our frequency analysis of HS 2201 +2610. We refer to the pulsation frequency identifications obtained by 
ISilvotti et alj ( l2002h and also list the mode degree (I) identification suggested there in Fig. 9. Calar Alto data: We estimate the 
errors in the frequency to be 0.2 yuHz (FWHM of the central peak in the window function) and having a detection limit of 1 . 1 mmi, 
which was calculated as four times the mean noise level in the residual DPT. We used (f4) and ( f5) in the prewhitening, because 
low-amplitude signals at these frequency ranges were also found in other runs (I Silvotti et alj|2002l) . OHP/LOI data: We calculate a 
frequency resolution of 0.54 //Hz (FWHM of the central peak in the window function) and a detection limit of 0.97 mmi (four times 
the mean noise level in the residual DPT). Prom the WPT Xcov23 runs, the OHP data from August 22/23 and the LOI data from 
August 29/3 1 and September 3/4 have been included. Mean sol ution: The frequenci es of the "mean" solution refer to BJD 2452788, 
which corresponds to the middle of the 7-year long data set in ISilvotti et alj ( l2007l) . Estimating the formal frequency resolution of 
the 7-year "mean" solution yields 0.0045 juHz. 
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mass main-sequence (M-dwarf) companion as the source of the 
54min signal can be excluded by geometrical considerations 
(Roche-Limits vs. separation in Keplerian orbits). The same is 
true for a sdB+WD configuration in a 54min orbit down to very 
low WD-masses (below 0.2 M0). In the case of deformation of 
the sdB by a hypothetical WD companion, causing ellipsoidal 
variation, the underlying orbital period would be 108 min, but 
we would then expect a signal at least ten times stronger down 
to low inclinations, dropping to the observed very low amplitude 
at only about 10°. A WD companion as the source of the signal is 
therefore also very unlikely. Por less massive and hence harder 
to detect companions it would be increasingly hard to prevent 
mass transfer in the system. Binary orbital motion as the cause 
of the long-period variation will immediately and conclusively 
be ruled out if its multi-periodicity is confirmed. 

The period and amplitude fall within the typical range for 
g-mode sdB pulsations, strongly suggesting that the variation is 
in fact due to a g-mode. Pinally, in comparison to the location 
of the other two known hybrid sdB pulsators in the log g - T^ff 
diagram of sdBV stars (see e.g. Pig. 1 in Lutz et al. 2008), the 
very similar parameters of HS 2201+2610 again imply that hy- 
brid behaviour may plausibly be expected. All three objects are 
located at the interface between the empirical /9-mode and g- 
mode instability regions. Hybrid pulsators have also been shown 
to exist in the partially overlapping regions of the /3 Cep / [SPB] 
stars , which resemble the p/ g-mode sdB pulsators in many ways 
(e.g. lHandler et al.'2004VJerzvkiewic z et al.l2005tlHandler et al.1 
I2OO6 ; Pigulski & Poimaiiski 2007a b). Hybrid behaviour is also 
observed in the less closely related 6 Set / yDor stars (e.g. 
iHenrv & Fekel 2005). 

From a theoretical point of view, modeling p-mod e sdB pul- 
sators has been very successful. Fontaine et al. (2006") report on 
the determination of the basic structural parameters of twelve 
p-mode sdB pulsators. Regarding the g-mode sdB pulsators 
th ere are still issue s . The "sec ond-generation models" presented 
bv lCharpinet etal.l (1 19971) and lFontaine et al.1 (^003) diverge by 
some 4000 K from the observations, i.e. they fall short of the 
empirical blue edge of the g-m ode pulsators. A solution to this 
problem has been suggested by iJefferv & Saiol (2006, 2001^: In 
their models, using Opacity Project (OP) instead of OPAL opac- 



ity tables and including Ni in addition to Fe to the opacity pro- 
file shifts the theoretical blue edge close to the empirical one. 
Both aspects of this opacity problem are, in close analogy, als o 
relevant for the /3Cep and [SPB] pulsators (iMiglio et al.ll2007h . 
In fact, the hybrid objects (both in the main-sequence as well 
as in the extreme horizontal branch pulsators groups) in partic- 
ular represent excellent test cases for the adequateness of the 
opacities used, and should provide strong constraints. Prom pre- 
vious experience, it appears to be clear that radiative levitation 
and self-consistent diffusion of the Pe/Ni group elements need 
to be included in the EHB models. It will be interesting to see 
if "third-generation models" incorporating both radiative levita- 
tion/diffusion and improved opacities will finally deliver an ad- 
equate description of the hybrid sdB pulsators while still repro- 
ducing the behaviour of the p-mode pulsators as consistently as 
before, and at the same time correctly describing the extent of 
the empirical g-mode instability region. 

5. Concluding outlook 

As recently published bv ISilvotti etall (l2007h . HS 2201+2610 
has a planetary companion which was detected using asteroseis- 
mological techniques. Long-term monitoring of the strongest p- 
mode pulsations made it possible to derive secular variations of 
the two main pulsation frequencies, leading to an evolutionary 
timescale of around 5-8-10* years, which matches theoretical 
evolutionary models of extreme horizontal branch stars. Periodic 
residuals in the O -C diagram (observed - calculated) analysis, 
attributable to light-travel time effects, revealed the presence of 
an msini - 3.2MJ^p companion around HS 2201+2610 in a 
3.2 yr orbit. We are in the process of continuing our regular long 
term photometric monitoring ofHS2201+2610 beyond the cur- 
rent coverage of two orbit cycles. In order to further characterise 
the companion, especially to determine its mass, it is necessary 
to constrain the orbital inclination of the system and to derive a 
more accurate asteroseimic mass of the host star. 

The mass of the host star is in principle accessible through 
asteroseismology, and Silvotti et al. (2002) already presented 
one possible solution. The determination of the fundamental pa- 
rameters of the star has however only weakly been constrained 
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in that work due to the small number of modes detected. Now 
that we have added HS 2201+2610 to the list of hybrid sdB 
pulsators, this finding might provide the additional constraints 
necessary for a reliable and accurate asteroseismic characteri- 
sation. At this point, we emphasize again the cuiTent need to 
improve on the model calculations first. For the individual ob- 
ject HS 2201-1-2610, it will be a special challenge to addition- 
ally drive ;?-mod es in an improved ve rsion of one of the possible 
models found bv lSilvottietan(l2002h . 
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